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Abstract 

Recent work has shown that non-local modifications of gravity involving terms such 
as irPRD~‘^R (and no cosmological constant) provide a phenomenologically viable al¬ 
ternative to ACDM. We first discuss the possibility that such non-local terms emerge 
in the far infrared from the running of a coupling constant associated to the R^ term 
in higher-derivative gravity, which, depending on the UV completion of the theory, 
can be asymptotically free in the ultraviolet and strongly coupled in the infrared. In 
this scenario the mass scale m of the non-local model emerges from dimensional trans¬ 
mutation, similarly to Aqcd for strong interactions, leading to a technically natural 
value and to a novel understanding of the scale associated to dark energy. Motivated 
by these findings, we then explore the possibility of generating strong infrared effects 
in Einstein gravity, with no Rf terms, as a consequence of the higher-derivative term 
generated by the conformal anomaly. 



1 Introduction 


At the fundamental level, quantum field theory is local. Still, in many situations non¬ 
locality emerges as a derived phenomenon. This can happen already at a purely classical 
level, when one integrates out some fast degree of freedom to obtain an effective theory 
for the slow degrees of freedom, or at the quantum level, when one considers the effective 
action that takes into account loop corrections involving light or massless particles. In 
recent years there has been a significant activity on effective non-local modifications of 
gravity, largely motivated by the aim of understanding the origin of dark energy. In 
particular, in [1] (elaborating on previous works related to the degravitation idea [2-4], 
as well as on attempts at writing massive gravity in non-local form [5,6]) we proposed a 
phenomenological modification of gravity, based on the non-local equation of motion 

- (l/3)m2 = SttGT^, . (1.1) 


The superscript T denotes the operation of taking the transverse part of a tensor (which is 
itself a nonlocal operation), □ is the covariant d’Alembertian computed with the curved- 
space metric and its inverse 0“^ is defined using the retarded Green’s function, to 
ensure causality. The factor 1/3 is a convenient normalization of the parameter m? in d = 3 
spatial dimensions. The extraction of the transverse part ensures that energy-momentum 
conservation is automatically satisfied. A closed form for the action corresponding to 
eq. (1.1) is not known. This model is however closely related to another non-local model, 
proposed in [7], and defined by the action 
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These two models are related by the fact that, linearizing over flat space the equations 
of motion derived from the action (1.2), one finds the same equations of motion as those 
obtained by linearizing eq. (1.1) over flat space. However, at the full non-linear level, or 
linearizing over a background different from Minkowski, the two theories are different. 

An intriguing aspect of the models (1.1) and (1.2) is that, in the comparison with cos¬ 
mological observations, they perform remarkably well. They both have the same number 
of parameters as ACDM, with the mass m replacing the cosmological constant. They dy¬ 
namically generate a dark energy and have a realistic background FRW evolution [1,7,8]. 
Their cosmological perturbations are well-behaved (a step that already ruled out several 
modified gravity models) and their quantitative effects are consistent with CMB, super¬ 
novae, BAG and structure formation data [9-11].^ This allowed a more detailed compar¬ 
ison with ACDM. Implementing the cosmological perturbations of the non-local models 
into a Boltzmann code and performing parameter estimation and a global fit to CMB, 
supernovae and BAO data, one finds that these models perform as well as ACDM, with 
comparable values of the [20] (in fact for model (1.1) the is even better than for 

^This should be contrasted with the non-local model proposed in [12-14]. This model does not involve a 
mass scale, and is rather constructed adding to the Einstein-Hilbert action a term of the form Rf{0~^R). 
The fnnction f(D~^R) is tnned so that, at the level of background evolution, this model closely mimics 
ACDM. One can then study its cosmological perturbations, and it has been found in [15] that the model is 
ruled out by the comparison with structure formation. Non-local long-distance modifications of GR have 
also been suggested in [16-19]. 
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ACDM, although not at a statistically significant level). Furthermore, parameter estima¬ 
tion provides a value of the Hubble constant Hq slightly higher than in ACDM, in better 
agreement with that obtained from the most recent local measurements [21], To the best 
of our knowledge, among the large variety of existing alternatives to the standard ACDM 
paradigm, the models defined by eqs. (1.1) or (1.2) are the only ones that are competitive 
with ACDM from the point of view of fitting current observations (at a level of accuracy 
which tests not only the background evolution but also the cosmological perturbations of 
the model), without being just an extension of ACDM with extra free parameters (see 
also [22]). Further conceptual and phenomenological aspects of these models have been 
discussed in [23-30]. 

These non-local models have been proposed on a purely phenomenological ground, 
and it is of course important to eventually understand how they could emerge from a 
fundamental theory. The purpose of this paper is to investigate the possibility that such 
terms emerge from infrared (IR) effects in gravity. We will first discuss the idea of gen¬ 
erating the dark-energy scale from the running of asymptotically-free coupling constants 
in E? extensions of Einstein gravity. We will then explore the possibility of generating 
similar infrared effects in Einstein gravity, with no E? terms, just as a consequence of the 
higher-derivative terms generated by the conformal anomaly. 

The plan of the paper is as follows. In Section 2 we recall known results on the 
running of the couplings associated to the 0{E?) terms. In sect. 3 we discuss the possible 
IR dynamics that could be generated by the running of these coupling constants, and the 
connection with the non-local models (1.1) and (1.2). Just as in QCD, the running of a 
coupling that becomes strong in the IR generates, by dimensional transmutation, a mass 
scale Arr, analogous to Aqcd for strong interaction. In sect. 4 we show that the dynamical 
emergence of this mass scale can have remarkable consequences on several aspects of the 
“cosmological constant” problem. An important role in the picture that we will develop 
is played by the conformal mode of the metric. Then, in sect. 5 we examine in more detail 
the dynamics of the conformal mode in E? extensions of GR. Einally, in sect. 6 we will 
consider Einstein gravity, without E? terms, and the possibility that similar effects arise 
in this case from the higher-derivative term provided by the anomaly-induced effective 
action. In sect. 7 we summarize our results. 


2 Non-local loop corrections in gravity 

2.1 Notations and conventions 

We consider the theory with action 


S 


*S'eh + Shd 

j 



{E — 2A) — (aiC*^ J- (I 2 E? J- , 


( 2 . 1 ) 


where mpi = l/(87rG)^/^ is the (reduced) Planck mass, Seh is the Einstein-Hilbert (EH) 
action with a cosmological constant, and S'hd is the higher-derivative term. Here 


~ + (1/3)R^ 
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is the square of the Weyl tensor, and 

rp _ d 2 


(2.3) 


-ARl, + R^ 


is the Gauss-Bonnet term. Observe that oi, 02 and 03 are dimensionless. We use the MTW 
sign conventions [31], so in particular = (—, -|-, -|-). The overall minus sign in front 

of the higher-derivative terms in eq. ( 2 . 1 ) is part of the definition of the a* coefficients, 
and is chosen so that these terms will appear with a plus sign in the euclidean action, 
see eq. (2.5) below. Loop computations are indeed typically performed with euclidean 
signature, and it will be important for us to be careful about the signs in the passage from 
Minkowskian to euclidean signature. The rotation to euclidean signature (-I-, -|-, -b) is 

obtained performing the Wick rotation, i.e. introducing euclidean time from = it. 
Then d^x —>■ —f((i^x)E, while —)• < 7 ^^^ and R transforms into the Ricci curvature 

computed with the euclidean metric, R^, without extra minus signs, and similarly for the 
Riemann and Ricci tensors. Then, 


(■S'eH + -S'HD)Mmk = -i {d‘^x)Eg^^'^ 


m 


PI 


(i?E ~ 2A) — (aiCg -|- (12Re T ®3-^e) 


• (2.4) 


Defining the euclidean action 5 euc 1 from = e ‘S'buci^ SeucI = we 

therefore find 


*S'euc1 = / {d'^x)E 


(2.5) 


We will henceforth drop the subscript E from euclidean quantities. Whether an action 
is written in Minkowskian or euclidean signature can be understood from the presence of 
the factor \/—g or respectively. We will generically denote by R^ gravity the theory 
containing the Einstein-Hilbert action plus all the terms R^, and E. 


2.2 gravity at low energies 

Let us first consider R^ gravity in an effective field theory approach, valid in the limit 
E <C Mpi. In this case the 0{R?) terms are simply seen as the first corrections that we 
meet as we approach the Planck scale. More precisely, if we linearize over flat space, g^y = 
dfiv + hfj,y, the Lagrangian density reads, schematically (i.e. displaying only the powers 
of hf^y and of derivatives, without writing explicitly any tensor structure nor numerical 
factors), 

£eh + >Chd ~ mpi(/r9^/i -b hdhdh + ...) — a [(nh)^ -b h{Oh)‘^ -b ...] ( 2 . 6 ) 


where, for the purpose of this schematic counting, a denotes generically oi or 02 (the 
Gauss-Bonnet term is a topological invariant and we will neglect it for most of the following 
discussion). In an effective field theory approach the degrees of freedom are read from the 
Einstein-Hilbert term, so in the spectrum we only have the massless graviton. Ganonical 
normalization is also performed with respect to the Einstein-Hilbert term, i.e. rescaling 
h^y —)• h^y/mp\, and therefore the Lagrangian density is rewritten as 




( hd'^h —^hdhdh + ... 

)- “ 

(□/i)2+ ^ h(nh)‘^ + ... 
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(2.7) 


3 













We see that the contribution to any process from the higher-derivative terms is suppressed, 
compared to the contributions coming from to the EH action, by a factor of order 
The standard lore is therefore that these corrections are irrelevant in the infrared, where 
E <C mpi- 

The idea that we want to explore in this paper is that, in an effective action that derives 
from a full quantum theory of gravity, ai and 02 are dimensionless coupling constants that 
run under renormalization group, so a* = ai{E). Depending on the sign of their beta 
functions, their absolute values can either grow or decrease as we run toward the IR. If a 
coupling grows in the IR, i.e. if the sign of its beta function corresponds to asymptotic 
freedom in the UV, we meet a situation conceptually similar to what happens in QCD. 
Running toward the IR this coupling will enter a strong coupling regime. The energy 
scale at which this happens defines a dynamically generated scale, analogous to Aqcd- 
At energies below this new scale, the functions ai{E) have a behavior which has nothing 
to do with the one computed in perturbation theory, and can in principle even develop 
singularities as E —?■ 0. In particular, if in this regime a 2 {E) oc E~‘^, in coordinate space 
we would have in the action a term Ra2{0)R oc RO~‘^R (see sect. 2.3 below), which would 
indeed reproduce the model (1.2). 

In this section we investigate this scenario. We will first discuss the running of the 
couplings associated to the 0{R?‘) terms, and we will then give arguments, drawing on the 
known behavior of singularities in the QCD Green’s functions, as well as on the dynamics 
of the conformal mode in R? gravity, that suggest that a IR behavior 02 ( 0 ) oc is indeed 
plausible. We begin by reviewing the known results on perturbative loop corrections to 
the ai couplings. 

2.3 Non-local form factors at one loop 

In general in QFT the running of coupling constants, which is more commonly expressed in 
momentum space, can also be expressed in terms of non-local form factors in the effective 
action. For instance, at the one-loop level the running of the electric charge in QED can 
be described in coordinate space by the effective action [32] (see also [33,34]) 

Se» = -)/ 

where 

1 1 

?{5) 

Here /U is the renormalization scale, e(^) is the renormalized charge at the scale /U and, 
for a single massless fermion, /3o = l/(127r^). The logarithm of the d’Alembertian can be 
defined for instance from 

=r 

Such non-local actions generate non-local effective equations of motion for the expectation 
values of the quantum fields. The corresponding equations of motion for the in-out matrix 
elements depend on the Feynman propagator and are acausal, but the effective classical 
equations of motion for the in-in expectation values, which can be computed using the 
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Schwinger-Keldysh formalism, involve the retarded propagator and are therefore non-local 
but causal [35,36]. In practice, using the in-in formalism turns out to be equivalent to 
performing the computation in euclidean space, and in the end rotating to Minkowskian 
signature while at the same time replacing the euclidean 0“^ with the Minkowskian 0“^ 
computed with the retarded Green’s function, as discussed in [32] and recently verified 
explicitly, for the case of E? gravity, in [37]. 

Non-local terms also appear in the renormalization of gravity with higher derivatives, 
as was first revealed using the heat-kernel technique in [32,38]. The result has also been 
checked with standard Feynman diagram computations [39,40]. These non-local terms 
appear when we consider gravity semiclassically and quantize massless scalars, massless 
spinors or massless vector fields over a fixed curved background. Furthermore, they also 
appear when we quantize gravity itself, as a result of the contribution of gravitons to 
the quantum loops. The contribution from matter field is very well understood, and 
summarized in textbooks such as [41,42]. In contrast, a consistent computation of the 
running due to graviton loops is a much more subtle issue. We examine these contributions 
separately. 


2.3.1 Loop corrections from matter fields 


To lowest perturbative order, the loop corrections induced by massless matter fields can 
be computed in semiclassical gravity, i.e. on a fixed curved background. Using 
and R^ as independent terms, the non-local correction to the euclidean action at one-loop 
takes the form^ 




1—loop 


NL 


2(47r)2 
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aR\og R + PRp^log 
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+jRp,p„log ] R^^^P^ 
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where □ is the generally covariant d’Alembertian in euclidean space. Switching 
basis of R‘^,C‘^,E, this expression can be rewritten as 


5 , 


1—loop _ 
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2 ( 47 r) 




( □ 


a Rlog ] R + P log ( — ) CP'^P'' 


+7 Ruupa log 


log 


RP'^ + R\og 
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( 2 . 11 ) 
to the 


( 2 . 12 ) 


where Cp^pa is the Weyl tensor, while the last structure corresponds to the Gauss-Bonnet 
term (2.3). The relation between the coefficients is [37] a = a+ ^+j, P = —2/3 — dy, 
7 = /3 -|- 7 . The value of the coefficients {q;,/3 ,7} can be found e.g. in [41] (since the 
logarithmic terms are fixed by the 1/e divergences of the effective action). In Table 1 we 
collect the values of {a, 7 } due to a massless scalar with a generic non-minimal coupling 

^ to the curvature, described by the euclidean action 


Ss = ^ / d'^xg^!'^ [gp'"dp^dy^ E iR(^^) , (2.13) 

^We follow the notation in [37], except that we define the constants extracting an overall factor 

l/[2(47r)^] in front of which makes easier the comparison with most of the literature, e.g. with 

refs. [39,40]. 
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a 


7 

scalar 

(l/2)p 

1/120 

-1/360 

fermion 

0 

6/120 

-11/360 

vector 

0 

12/120 

-62/360 


Table 1: The value of the parameters {a,/3,7} due to loops of massless matter helds. 


as well as the contributions from massless spinors and massless vector fields. The value of 
the scalar held is given in terms of the parameter ^ = ^ — (1/6), and = 0 corresponds to 
a conformally coupled scalar held. 

The cosmological effect of these non-local terms has been recently studied in [37]. 
Even if enhanced by a logarithmic factor, a term R\og{—U/^’^)R can compete with the 
Einstein-Hilbert term mpjii only for nearly Planckian curvatures. Therefore, the inclusion 
of the UV form of the loop corrections can be relevant for addressing issues such as the 
resolution of cosmological singularities near the big-bang, or the emergence of the classical 
behavior from the full quantum theory, which are the issues addressed in [37].^ Here we 
are rather interested in the opposite limit of curvatures much smaller than mpi and very 
large distances, i.e. in the far IR, with the aim of understanding the origin of dark energy. 

It is also instructive to compute the one-loop renormalization of the R?‘ and terms 
induced by massive (rather than massless) matter fields in a curved background, to un¬ 
derstand the interplay between non-locality and the mass of the field. This computation 
has been performed in [39,40] for scalar, spinor and vector fields (neglecting the renormal¬ 
ization of the Gauss-Bonnet term), using a mass-dependent renormalization scheme that 
allows one to correctly recover the decoupling of massive helds in the limit of small mo¬ 
menta. Consider for instance a real scalar held with mass m<j and a generic non-minimal 
coupling described by the euclidean action 

= d^,(t)dy(t) + ml4>^ + iRcj)^) . (2.14) 


The corresponding one-loop contribution to the euclidean effective action is [39,40]^ 


^1—loop 
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1 + kR{n) 
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(2.15) 


^See also [43], where one-loop corrections to the theory computed expanding over de Sitter space, 
and involving corrections of the form log R/, are used to construct an early-universe inflationary 
model. Non-local extension of the R^ theory relevant in the UV have also been studied, with different 
motivations, in several papers, see e.g. [44,45]. 

'^This result was hrst obtained in [39] using a non-covariant computation based on an expansion of the 
action over flat space, and it was then checked in the same paper using the covariant heat-kernel technique. 
In the covariant technique the quantity which is computed is dehned by def 

expfr^'^^}. Since the euclidean action enters as we have This is the origin of the 

overall sign difference between the result of the covariant computation (eqs. (2.15) and (2.18) of [39]) and 
of the non-covariant computation, eq. (3.9) of ref. [39]. I thank E. Gorbar for clarifying me this point. 
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where 1/e = 2/(4 — D) + log{4:TTfj?/m'^) — 7e comes from the dimensional regularization 
scheme in D space-time dimensions, and the form factors kwi^) and are given by 


where 


kw{^) 

kR{a) 


8^ 
15 a4 

1 




108 o2 


A = 1 - - 
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2 1 
45 a2 IM ’ 
2A A 1 \ 
3a2 “ ISy 
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~ 2160 ’ 




1 1 \ 
1 ^^ 144 ; 


4n 

□ — 4m2 


(2.16) 

(2.17) 

(2.18) 


A few comments on this result are in order. First, observe that the first two terms in 
eq. (2.15) correspond to the renormalization of the cosmological constant and of Newton 
constant, respectively. Even if they depend on the renormalization scale [i (through the 
parameter e), they have no dependence on the □ operator and, in this sense, no running in 
momentum space. This is a trivial consequence of the fact that DA = 0, while Di? is a total 
derivative, so we cannot obtain a form factor by acting with the □ operator on A or on R. 
In contrast, the and the R? terms acquire corrections of the form Cpupakw{^)C^'^^^ 

and RkR{0)R, in full analogy with eq. (2.8). In the UV limit D/m^ ^ 1 we have (see 
also [33,46]) 


kw{^) — 

kR{n) ~ 
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23 


--log^ + —+ 0 ^log 


60 
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mt 
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-F^log-^+ ( - 
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mt 
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+ ^ + 

1080 18 


f)- 

+0 



(2.19) 

( 2 . 20 ) 


Therefore, after subtracting the divergent parts, in the limit D/m^ 1 we remain with 


( 2 . 21 ) 

(plus local terms independent of □) in full agreement with eq. (2.12) and Table 1 (the 
renormalization of the Gauss-Bonnet term was not included in this computation). Observe 
also that, in the conformal case mg = 0,.^ = 0, the form factor A;_R(n) becomes local, 
A;R(n) = —1/1080, in agreement with the prediction that can be obtained directly from 
the conformal anomaly.^ 

®Again, checking the sign is a bit tricky due to the many different conventions in the literature. The 
standard anomaly result can be read from eq. 6.124 and Table 1 on page 179 of the Birrell and Davies 
textbook [41]. One must however take care that Birrell and Davies use the signature (-1-,—,—,—) and 
their definition of Riemann tensor is opposite to MTW. To go to MTW sign conventions we therefore must 
switch gfj,u —> —Qiiu (so □ —^ —□) and R — > -\-R. Taking into account that this flips also the sign of the 
trace , we find that the correct result for the anomaly of a single massless conformally coupled scalar 
field, in MTW conventions, is T = + i 8 o( 47 r^) From the signature (—,-I-,-b,-b) we can then trivially 
rotate to euclidean space without getting extra minus signs. In eq. (3.12) of [39] the sign of the anomaly 
in the euclidean space was incorrectly taken to be a minus, and the result was checked using instead 
see footnote 4, so these two sign errors compensated. 


^1—loop _ 
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In contrast, in the opposite limit D/m^ <C 1, eqs. (2.16) and (2.17) give 


kw{a),kR{a) = 0{a/m^), ( 2 . 22 ) 

corresponding to the decoupling of particles with mass large compared to the momentum 
scale, which is explicit in the mass-dependent subtraction scheme used in [39,40] (see also 
the discussion in [47]). Furthermore, beside being small, a term 0(n/m^) is local. This 
shows explicitly that non-local terms appear in the effective action as a consequence of 
the running in the loops of particles that are either massless, or light with respect to the 
energy scale considered. 

2.3.2 Gravitational loop corrections in Einstein gravity 

In semiclassical gravity, oi, 02 , 03 are simply parameters of the effective action. When we 
go beyond the semiclassical approximation, and treat gravity itself at the quantum level, 
ai and 02 become genuine coupling constants, and the perturbative expansion is organized 
in terms of the three coupling of the theory G, ai and 02 (while the Gauss-Bonnet term 
is a topological invariant, and will be neglected in the following). Their renormalization 
group equations also receive contributions from the purely gravitational sector, i.e. from 
graviton loops. However, these contributions depends crucially on the UV completion of 
the theory, as we now review. 

If one starts from the Einstein-Hilbert action without higher-derivative terms and 
computes the one-loop corrections, one in principle finds and terms, as was first 
shown in the classic paper by’t Hooft and Veltman [48]. The original computation of [48] 
gives, in our notation, a = 90/360 = 1/4 and /3 = 126/360 = 7/20, which are the values 
quoted in [37]. However, these divergences have no physical meaning, since they can be set 
to zero using the equations of motion [48], so that pure gravity at one loop is finite on-shell. 
Equivalently, they can be set to zero with a field redefinition —)■ g^u+(-{ciR^iv+C 2 Rgij,v)- 
More generally, if one quantizes the pure Einstein-Hilbert action, the coefficients of the 
R^ and R^^^, terms which are generated at one loop depend on the gauge used, and one 
can find gauges in which these divergences are absent, so that the theory is one-loop finite 
even off-shell [49]. Thus, in pure Einstein-Hilbert gravity there is no sense in which we can 
compute the gravitational contribution to the running of the coupling constants associated 
to the R^ terms. The same is true if we consider the higher derivative theory (2.1) in an 
effective action approach, since the higher-derivative term add extra vertices, but does not 
change the structure of the propagator. 

2.3.3 Gravitational loop corrections in Stelle theory 

The situation changes in an interesting manner if, rather than using an effective action 
approach, in which the higher-derivative term is seen as a correction valid only in the low- 
energy limit, we try to take the action (2.1) as a full UV complete theory, and quantize 
it. In this approach, the higher-derivative theory (2.1) is usually known as Stelle theory. 
In this case the propagator is read from the full quadratic term in eq. (2.6) and, again 
schematically, is of the form 

1 

rrip^k'^ + ' 
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G{k) 




(2.23) 



(before any canonical normalization) so that now it goes as 1/k^ in the far UV. This 
improved behavior of the propagator allowed Stelle to prove that this theory is renormal- 
izable [50].® Having a renormalizable qnantum theory of gravity is quite remarkable. This 
however comes at a very high price. Namely, the spectrum of B? gravity, which is now 
read from the full propagator, consists of the usual massless graviton, a massive spin-2 
particle and a massless spin-0 particle [52,53]. The massive spin-2 state has negative 
kinetic energy, so is a ghost, and its squared mass is 


9 m 

1712 = 


PI 


2ai 


(2.24) 


The massive spin-0 state has a normal kinetic term, and squared mass 

( 2 . 

and is a tachyon if 02 > 0. Because of the ghost, and (in case) of the tachyon, 
gravity, treated non-perturbatively as a full UV complete theory, is not a consistent theory 
(otherwise, the problem of finding a consistent theory of quantum gravity would have been 
solved in the late 1970s...). There have been attempts at exorcising the ghost based on 
the idea that radiative corrections shift its pole into the complex plane [54,55] , so that the 
ghost becomes unstable and does not appear in the asymptotic states. To leading order 
in 1/V, where N is the number of conformally coupled matter fields, the theory turns 
out indeed to be unitary. However, attempts at proving unitarity beyond leading order 
based on the gauge-dependence of the position of the ghost poles fail [56], so what happens 
to all orders is basically unknown. In any case, there is a long tradition of ignoring the 
ghost problem and see what the R? theory has to say at a full non-perturbative level, in 
the hope that this will at least help to unveil properties of a fully consistent theory of 
quantum gravity. In particular, one can compute the running of the coupling constants 
associated to the 0{R?) terms. As we mentioned above, in Einstein-Hilbert gravity this 
computation leads to beta functions that are not even gauge-independent. In contrast, in 
Stelle theory the computation is well-defined, and the beta functions associated to the R?, 
and Gauss-Bonnet terms are physical and gauge-independent (while, for the Newton 
constant and the cosmological constant, only the dimensionless combination GA. has a 
gauge-independent beta function) [53,57,58]. For Stelle theory, the computation of the 
gravitational contribution to the beta functions of the at coupling was first performed 
correctly in [58] (building on earlier work in [53,57]), and the result was confirmed in [59]. 

®Of course, the full story is more complex. The B? and terms also lead to vertices growing as k'^, 
rather than kF' as in Einstein-Hilbert gravity. In GR, since the propagator go as 1/kF and the vertices as 
k^, the superficial degree of divergence is D = 4L — 21 + 2V, where L is the number of loops, I of internal 
lines, and V of vertices. Combined with the topological relation L = 1 + I — V this gives D = 2L + 2, 
so the degree of divergence grows with the number of loops. In contrast, for Ehd, D = 4L — 4/ -|- 4T 
which, together with L — I + I — V, gives D = 4, independent of the number of loops [50]. This improved 
behavior is at the basis of Stelle’s proof of renormalizability. Observe also that, for renormalizability, it 
is crucial to include both the and the terms in the action, since they contribute separately to the 
propagator of the spin-0 and massive spin-2 excitation of the theory. Otherwise, the propagator of one of 
these excitations would still go as 1/fc^, while the vertices grow as k^, leading to an even worse UV behavior 
than in Einstein-Hilbert gravity. Furthermore, even the Gauss-Bonnet term is required for multiplicative 
renormalizability (see the recent discussion in [51]). 


mg = - 
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The result is conveniently expressed trading ai and 02 for two couplings A and u defined 
by 

so that the Minkowskian action (2.1) reads (neglecting the Gauss-Bonnet term and the 
cosmological constant) 


S = 



m 


PI 




(2.27) 


The corresponding beta functions are given in eqs. (19)-(24) of [58]. In particular, in¬ 
troducing the logarithm of the energy-scale parameter t = (dvr)”^ logi?//rO) the purely 
gravitational contribution is 

dX-^ 
dt 
du 
dt 

where /3g = 133/10 is positive. It is important to stress that these beta functions have 
been computed in the UV limit E ^ mp\. In the regime E <C mpi which is interesting 
for us, these equations will receive corrections from the Einstein-Hilbert term, which have 
not been computed. The solution of eq. (2.28) is X~^{t) = X~^{t = 0) -|- fdgt, i.e. 


= A 


g > 


> /lO 2 183 

= -My” +io“+ 
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(2.28) 

(2.29) 


l + A(^)/3,log(E/^)- 

where fdg = ;0g/(47r)^. Whether this corresponds to a function X{E) that becomes large 
in the UV and small in the IR, or viceversa, depends of course of the sign of A(/r). Since 
Pg > 0, if A(/x) > 0 we have X{E) —>■ 0 as E —)• 00, and conversely X{E) gets large and 
formally diverges at a hnite energy scale in the IR. The opposite situation takes place if 
A(//) < 0. Observe from Table 1 that the contribution to the beta function of ai = X~^ 
coming from each types of matter fields has the same sign as the gravitational contribution, 
so the same behavior takes place in the full theory with gravitational plus matter loops. 

The evolution of u is instead determined, at least in the super-Planckian regime in 
which eq. (2.29) holds, by the two roots of the equation (10/3)u^-|-(183/10)u-|-(5/12) = 0, 
which are given by ui ~ —5.46714, U 2 — —0.02286. Writing eq. (2.29) as 


du 10 , . 

— =- — X{u-ui){u-u 2 ), (2.31) 

we can easily read the sign of du/dt and the direction of the RG flow, for different values 
of u, and a given sign of A. For A > 0 the one-dimensional flow in the u variable is 

depicted in the upper panel of Fig. 1 and for A < 0 is shown in the lower panel. Observe 

in particular that, for A > 0, u = tt2 is a UV fixed point, which attracts all initial values 

u{t = 0) > ui. For A < 0 the value u = ui is a UV hxed point at its attraction basin is 
given by u{t = 0) < U 2 - 

One might wonder how is it possible that a computation which is well-defined in the 
full B? theory is not well-defined in Einstein-Hilbert gravity, since the latter is formally 
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X>0 
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X<0 


——•-— 

Ui U 2 0 u 

Figure 1: The RG flow of u{t) for A > 0 (upper panel) and A < 0 (lower panel). The 
direction of the arrows indicate the flow of n as t increases, i.e. as we move from the IR 
to the UV. For a given initial value in the UV, the flow as we run toward the IR regime 
is therefore obtained reversing the sign of the arrows. 
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Figure 2: The graviton propagator and a one-loop correction. 


obtained taking the limit ai, 02 —)• 0 in the action (2.5). The reason is that, when we take 
the theory as a full non-perturbative theory, the structure of the divergences, which 
determine the beta functions, changes discontinuously if we set oi or 02 to zero. Indeed, 
consider for instance a loop correction to the graviton propagator such as that shown in 
Fig. 2. Even when the momentum p of the external line is sub-Planckian, \p‘^\ <C rripj, 
in the loops we will have propagators snch as G{k), given in eq. (2.23), and G{k — p), 
where k is the loop integration variable and p the external momentum. The divergence is 
determined by the behavior of the propagator (and of the vertices) for A: —>■ 00 . For any 
a 7 ^ 0 the propagator at A: —?• 00 behaves as l/(aA^), even when the external momenta 
are sub-Planckian, while for a = 0 it behaves as 1/A^ and the structure of the divergences 
is completely different. In other words, the limit a —)• 0 does not commute with the 
integrations over d'^k in the loop integrals. So, despite the fact that 5 hd is a higher- 
derivative term, which naively should be irrelevant at low energies, the structnre of the 
divergences of Stelle theory is completely different from that of GR, even when the external 
lines are sub-Planckian, and the limit a —)• 0 (or, more precisely ai —?■ 0 or 02 —>• 0) of 
the beta functions of gravity is not smooth. This fact has already been stressed in 
sect. 3.1 of [53], where in particular it was observed that the beta function of the pure 
Weyl theory cannot be obtained setting 02 = 0 in the beta function of the theory with a 
higher-derivative term aiC^ -|- a 2 R^- 
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2.3.4 Gravitational running in a consistent theory of quantum gravity 

The result (2.28, 2.29) cannot be applied directly to our problem, namely the running of 
02 at low energies in the actual universe, for two reasons. First, the computation leading 
to eqs. (2.28) and (2.29) has been performed in the deep UV region, E ^ mpi. In the 
regime E <C mpi, which is interesting for us, these equations will receive contributions 
from the Einstein-Hilbert term, which have not been computed. This is just a technical 
problem. However, a second and more fundamental problem is that, barring a resolution 
of the ghost problem, Stelle theory can anyhow only be taken as a toy model, rather than 
a theory describing Nature. 

However what we learn from this toy model is that, in a theory of gravity with a 
good UV behavior, the computation of the beta functions associated to the Oj couplings 
is well defined, and these quantities run with energy. In the end we expect that, whatever 
will be the consistent fundamental quantum theory of gravity, in that context it will be 
possible to compute consistently the running of the coupling constant associated to the 
and term in the corresponding low-energy effective action, including the gravitational 
contribution. ' 

From the above analysis we also learn that the running of the couplings ai{E) depends 
on the precise form of the UV completion of the theory, even at E <C mpp Indeed, 
the running of ai[E) computed in a consistent UV completion of GR cannot reduce at 
low energies to that computed in the effective low-energy theory (i.e. in the Einstein- 
Hilbert action supplemented by the higher-derivative terms, considered as corrections in 
the effective held theory sense), simply because there is no consistent result in the low- 
energy effective theory. The beta functions computed in the low-energy effective action 
based on the Einstein-Hilbert term, in which the R? terms are treated as corrections, are 
not even gauge-invariant. This of course does not happens e.g. in QCD, because QCD is 
renormalizable, which in the end means that, even if we eventually embed it in a larger 
theory, the details of the UV completion, much as the details of the regularization scheme, 
are irrelevant at low energies. This is not the case for Einstein-Hilbert gravity, which is not 
renormalizable. In a sense, this phenomenon is an interesting example of UV/IR mixing, 
since it implies that the IR physics of GR can sensitive to the UV completion. 

Thus, the points that we will take home from the following analysis are the following. 

1. In a fundamental and consistent quantum theory of gravity the couplings associated 
to the R? and term will in general run with energy, and their beta functions will 

^For instance, in string theory one could in principle compute the string-loop correction to gravitational 
scattering amplitudes directly at the world-sheet level, and reconstruct from it the corresponding loop 
corrections in the target-space effective action. Actually, in the effective action of type II string theory 
in 10 dimension there is no E? term because of supersymmetry. However, in general, the presence and 
the structure of 0(B?) corrections that one obtains in four dimensions will depend on the type of string 
theory considered, as well as on how supersymmetry is broken in the compactihcation to four dimensions. 
Observe that it is sometimes argued that, since the fundamental string theory is ghost-free, in the low- 
energy effective action can appear only the ghost-free combination given by the Gauss-Bonnet term. This 
argument is however wrong. Being ghost-free is a property of the full theory, and if we write the effective 
action derived from a ghost-free theory in an expansion in derivatives and truncate this expansion to a 
finite order (e.g. in our case keeping only the terms up to 0{R^)) apparent ghost state will in general 
appear, with a mass of the order of the cutoff scale of the effective theory (in our case mpi or the string 
mass scale), but they are just an artefact of the truncation. See e.g. [60] for a nice explicit example. 
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well-defined and gauge-independent.® 

2. The precise form of the corresponding RG equation, even in the IR regime E <C mpi, 
depends on the specific UV completion of the theory, and is therefore at present 
completely unknown. We cannot use the results from GR plus higher-derivative 
terms, considered in an effective action approach, because the corresponding beta 
functions are non even gauge-invariant. We cannot use the results of Stelle theory 
(2.28) and (2.29), either, because, first of all, these results only hold in the regime 
E ^ mpi and, more importantly, Stelle theory anyhow is not a consistent UV 
completion of gravity because of the ghost. 

R is also interesting to observe that, to obtain the model (1.2) in the IR, it is not 
a priori necessary to match the sign of 02 (r) in the UV with its (positive) sign in the 
IR. Indeed, Fig. 1 gives an explicit example of a RG flow that changes the sign of u, 
and therefore of 02 . A renormalization group flow that changes the sign of the coupling 
02 ( 0 ), flowing from a value 02 (^) < 0 at some scale /r corresponding to the energy scale 
of inflation, to a value 02 (^ 1 ) = +A^j^/n^ in the IR (where App is a constant that we 
will discuss below), would be particularly intriguing since it would interpolate between 
the Starobinsky model [61] in the UV and the non-local model (1.2) in the IR, thereby 
providing a unified description of inflation in the early universe and dark energy in the 
present epoch.® 


3 Infrared dynamics 


Even if we do not know the form of the IR running of the couplings ai{E) in the funda¬ 
mental theory of gravity, we can still analyze some general scenarios. In the end, there 
are basically two natural possibilities. Either |a 2 (E')| grows in the IR, or it decreases. The 
interesting situation, from our point of view, takes place when in the IR |o 2 (E)| grows 
and, at least at sufficiently low energies, 02 (E) > 0 (since this sign allows us to eventually 
match 02 ( 0 ) with the sign of the model 1.2). In this case, the situation will be quite 
similar to that in QCD. To stress the similarity, let us use the notation a 2 {E) = g^{E) 
and assume that g^{E) is asymptotically free in the UV so, as long as we are in the 
perturbative regime. 


with oq > 0. This implies that g^{E) eventually becomes large in the IR. Similarly, we 
can write ai(E) = f‘^{E) and study the scenario in which f‘^{E) is asymptotically free in 
the UV, and 


f^(E\ = /^(r) 

1 -h/2(/i)/3olog(E/^) ’ 


(3.2) 


with (5q > 0. To understand what happens when we enter more deeply in the IR regime is 
of course a difficult problem. As already observed recently [51,62] a pure B? theory, with 


®Of course, in some special case, the beta functions can be zero, e.g. if they are protected by super- 
symmetry. 

thank Eugenio Bianchi for a discussion from which this idea emerged. 
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no Einstein-Hilbert term, has a confining Newtonian potential V{r) oc r, because of the 
l/|k|^ behavior of the propagator in the static limit. Its infrared limit is therefore very 
different from Einstein gravity, and indeed one can even expect that in this theory gravitons 
should be confined. When we also add the Einstein-Hilbert term, at energies E <C m-p\ 
the standard term from Einstein-Hilbert action dominates over the term and 

we recover the standard behavior of gravity. Still, in the infrared, from the point of view 
of the quantum theory of gravity, we are in a peculiar situation: loops associated to an 
expansion in powers of Newton’s constant are suppressed by powers of GE'^ ~ <C 1 

and are therefore irrelevant. However, if the the couplings ai{E) and/or 02 (E) run toward 
a strong coupling regime and formally diverge at some energy in the IR, quantum loops 
associated to them can be important. As a result, the form factors associated to the 
and terms can become quite different in the IR. To determine the expression of 
these form factors in the IR is of course a difficult non-perturbative problem. However, 
experience with QCD allows us to understand some aspects of it. In particular, we know 
that the logarithmic running of an asymptotically free coupling constant generates by 
dimensional transmutation a renormalization group invariant mass scale, that we denote 
by Arr, analogous to Aqcd-^° More precisely, once we fix observationally the value of 
at some scale /i, formally, as we run toward the IR, the one-loop expression for 
/^(E) will hit infinity at some scale Aj given by 

Similarly, the one-loop expression for g^{E) will hit infinity at some scale 

Ag = //exp|-• (3-4) 

As we run toward the IR, the first scale that we hit, where one of the two running couplings 
formally diverges, defines the strong-coupling scale of the theory, so Arr = max(Aj, A^).^^ 
We therefore have a mass scale that appears in the IR regime by dimensional transmu¬ 
tation. Then, in the strong coupling region, = /^(D/App) and = g^(n/A^j^). In 
principle, it is possible to obtain a power-like dependence of and on □, from a 
resummation of leading logs, of the form^^ 


00 


E 

n=0 


n\ 



n 



(3.5) 


similar scenario has been discussed recently, in [51,63], although in these cases dimensional trans¬ 
mutation is rather used to generate the Planck scale from a pure higher-derivative gravity. The idea has 
been around in different forms since a long time, see e.g. [64,65]. In ref. [66] similar ideas have been applied 
to the coefficient of the Gauss-Bonnet term, using arguments from supergravity. 

^'^One might wonder whether one of the couplings could still remain weakly coupled in the regime where 
the other is strongly coupled. The answer is in general no. For instance, a pure Weyl-square theory, which 
only has the coupling, is not multiplicatively renormalizable, and generates a term already at the 
two-loop level [67]. Therefore, when enters the strong-coupling regime, it necessarily induces a strong 
coupling in g^. For this reason, even if for instance the coupling 02 associated to the term should not 
have the appropriate sign for becoming strongly coupled according to its own one-loop RG equation, it 
could eventually be driven toward strong coupling by the coupling ai. 

^^See [68] for similar considerations applied to the running of Newton’s constant, and ]69,70] for earlier 
related ideas on how the large-distance behavior of quantum gravity might affect the cosmological constant. 
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for some value of v. The model (1.2) would then be obtained if = 2. Again, to determine 
the precise form of g^{^) in the IR is of course a difficult perturbative problem. However, 
it is interesting to observe that in QCD the introduction in the action of a non-local term 


m 




1 


-| ab 


fiiy 5 


(3.6) 


(where is the covariant derivative and m is a mass scale) correctly 

reproduces the results on the non-perturbative gluon propagator in the IR, obtained from 
operator product expansions and lattice QCD [71-75] (see also the review [76]). In the 
gauge = 0 we have 


m 

T 



■ 1 ■ 


[d'P 


ab 




+0(A3) 


(3.7) 


and therefore (3.6) is a mass term for the gluon. Observe that the expression (3.6) is non¬ 
local but gauge-invariant. Physically, the introduction of this term in the effective action 
can be seen as describing the emergence of a dimension-2 gluon condensate, (A“A“) ^ 0. 
The same operator has been proposed in 3-dimensional YM theory [77,78], to describe 
the generation of a mass gap due to IR fluctuations. 

The situation is quite similar for the term Rn~‘^R in the gravitational case. Indeed, 
consider the conformal mode (t{x) of the metric, defined choosing a fixed fiducial metric 
and writing 

g^,Ax) = . (3.8) 

Using for simplicity a flat fiducial metric the Ricci scalar computed from the 

metric g^y = is 

R = —66“^°^ (□u -|- dfj,ad^a) . (3.9) 

Therefore, to linear order in cj, 

R=-Qna + 0{P), (3.10) 

and (upon integration by parts) 

r\r = ^qP + 0{P). (3.11) 


The P term is just a mass term for the conformal mode, while the higher-order interaction 
terms on the right-hand side of eq. (3.11) (which are non-local even in a) are required 
to reconstruct a diff-invariant quantity.In this sense the term Rn~^R is completely 
analogous to the term (3.6) in the Yang-Mills case, and its introduction describes a sort 
of “condensation” of the conformal mode. The fact that the term RO~'^R has this special 

Actually, there is here a subtle technical point, both in the gravitational and YM case. Indeed, in 
eq. (3.11) we have implicitly using both the equality = 1 , and = 1. While the former is 

correct by definition of the latter is not true since □ has a non-trivial kernel, given by the fnnctions 
/ such that □/ = 0. On such functions the formal equality = / is not true. Thus, to be accurate, 

we should rather write in eq. (3.11) aO~^Oa instead of cr^. Note that the expression aD~^Oa has an 
invariance under constant shifts of the conformal factor, cr —>■ cr + c, which is lost if we write it as a^. The 
equality aO~^Da — becomes however correct for all functions a such that Da 7 ^ 0. 
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physical meaning, and the analogy with the YM case, makes it more plausible its emergence 
in the IR limit of E? gravity. 

Observe also that the requirement of general covariance does not fix uniquely the term 
that reduces to to quadratic order in a. Indeed, as we already mentioned, the equations 
of motion of the the models (1.1) and (1.2) have the same expansion to linear order over 
Minkowski space, so the corresponding actions are the same in an expansion to quadratic 
order over Minkowski, and they both reproduce the term oc cr^ in the action, although with 
different non-linear completions. Both models can therefore be in principle justified on 
the basis of the argument that a mass should be dynamically generated for the conformal 
mode. 

These arguments therefore suggests that, in the IR 

a2(n)~^, (3.12) 

where Arr is a dynamically generated mass scale, analogous of the QCD Lambda param¬ 
eter Aqcd- The parameter Arr is related to the mass scale m that appears in eq. (1.2) 
by 

^RR = ■ (3.13) 

As for the numerical value of Arr, just as in ACDM the energy fraction Ha is a derived 
parameter, fixed by the flatness condition, in the non-local model the parameter m is 
similarly fixed, again by the flatness condition. In ACDM one finds in this way that, to re¬ 
produce the observations, the cosmological constant must be of order Hq] not surprisingly, 
in the non-local model one finds that m? Hq [1,7,8]. Thus, from eq. (3.13), 

A^r = O(Fompi). (3.14) 

More precisely, using the best-fit values m ~ 0.67L7o for model (1.1) [1], or m ~ 0.28iLo 
for model (1.2) [7] we get Arr ~ 0.4(i7o”^Pi)^'^^ and Arr ~ 0.3(L7o”ipi)^'^^; respectively. 
Thus, numerically, 

ARR~0.8meV [model (1.1)] , (3.15) 

and 

ARR~0.5meV [model (1.2)] . (3.16) 

Finally, we should also in general take into account the effect of the running of the Weyl- 
squared term in the IR. The argument on the dynamical mass generation for the conformal 
mode says nothing about the behavior of the form factor for the Weyl-squared term, since 
the Weyl tensor vanishes for a conformally fiat metric. If the form factor /^(Fl) associated 
to the Weyl-squared term has no significant IR enhancement, and more generally if it 
grows slower than l/D^ in the far IR, its contribution to the cosmological dynamics will 
be negligible compared to the Rn~‘^R term, and we end up exactly with the non-local 

^"^However, in B? gravity the interpretation in terms of the dynamics of the conformal mode is compli¬ 
cated by the fact that the conformal already gets a kinetic term, as well as a mass term, at tree level from 

the B? term, as we discuss in sect. 5. This will eventually lead us to consider, in sect. 6, the dynamics of 
the conformal mode in Einstein gravity supplemented by the anomaly-induced effective action, where no 
such complication arises, and one can really have a dynamical mass generation for an otherwise massless 
mode. 
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models of the form (1.1) or (1.2) (depending on the precise form of the non-linear terms, 
since the two models are equivalent at linear order over flat space). In contrast, if also 
/^(□) develops a singularity proportional to a corresponding term proportional to 
should be included in the cosmological model. This would not affect the 
background FRW evolution, since vanishes in FRW, but would give a contribution 

at the level of cosmological perturbations. 


4 A solution to the naturalness problem 


The above mechanism for the emergence of a term rin?R^~^R in the long-distance effective 
action gives a new perspective on the naturalness problem of the cosmological constant. 
As is well known, there are two aspects in the cosmological constant problem: the “old” 
cosmological constant problem, namely why we do not observe a huge value for the vac¬ 
uum energy; and the “new” cosmological constant problem; namely, assuming that some 
mechanism sets to zero the vacuum energy, what fixes in a technically natural way the 
observed energy scale associated to dark energy? Here we have nothing to add to the 
“old” cosmological constant problem, since we have simply fixed the renormalized value 
of the cosmological constant to zero.^^ In contrast, concerning the naturalness problem 
of the scale associated to dark energy, the above scenario provides a new and interesting 
viewpoint. In the mechanism that we have proposed, in R?‘ gravity the scale Arr emerges 
from dimensional transmutation, similarly to Aqcd- There is no issue of naturalness for 
such a quantity, which is determined by the logarithmic running of a dimensionless cou¬ 
pling constant. Of course, the value of the scale generated dynamically in this way cannot 
be predicted, just as we cannot predict the value of Aqcd- The fact that Arr turns out 
to be of the order of a milli-eV is simply a fact of nature, just as Aqcd turns out to be 
of order 200 MeV. However, quantities such as Aqcd and Arr do not have a power-like 
dependence on the masses of the particles in the theory, and are therefore technically 
natural. 

Another attractive aspect of the above construction is that, to explain dark energy, we 
do not need to introduce a fundamental particle with an astonishingly small mass m ~ 
Hq ~ 10“^^ eV, as is the case for instance in massive gravity. In the above construction 
the fundamental mass scale of the theory, emerging from dimensional transmutation, is 
Arr, and in terms of it the effective action in the IR reads 


^Mink — 






(4.1) 


Even if, just as in QCD, the value of Arr cannot be predicted by the model, still the 
numerical value required by the comparison with the observation, which turns out to be 
of the order of the meV, does not look particularly surprising. An extremely small mass 
scale only emerges when, in eq. (4.1), we factorize mpj/2, rewriting the action as 


•^Mink — 




^^RR d ^ p 

2 A^2-^ 


m 


PI 


(4.2) 


mechanism that could naturally set to zero the value of the vacuum energy density is proposed 

in [79]. 
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(4.3) 


so the mass scale m defined in eq. (1.2) is given by 

mpi 

which is suppressed with respect to Arr by an extra factor of order ApR/mpi, and is then 
of order 10“^^ eV. However, the fundamental mass scale generated by the theory is App, 
not m. 


5 Conformal mode dynamics 

It is interesting to understand in more detail the consequences on the conformal mode of a 
RG flow that runs from aiB? at high energies, with oi = aiip) a constant, to 
in the IR. We expands on the technical detail in app. A, and here we summarize the main 
results. First of all, let us linearize for simplicity over flat space, and restrict to the scalar 
sector of perturbations. Then the metric can be written in terms of the gauge-invariant 
Bardeen potentials and T as 

=-(l-2^)dt2 + (i + 2$)dx^ (5.1) 

Equivalently, we can introduce the conformal factor a and another independent perturba¬ 
tion (fi, defined writing the metric perturbation in the scalar sector in the form 

ds^ = [—(1 — 2ip)dt^ + (1 + 2(/?)dx^] . (5.2) 

At the level of linearized theory, = 1 + 2a + 0{a‘^), so 2it = <I> — T and 2f/? = <I> -|- T. 
In GR, without terms, of course only the tensor modes are dynamical, and scalar 
perturbations are non-dynamical constrained fields that satisfy a Poisson equation, rather 
than a Klein-Gordon equation. Indeed, in the presence of matter with energy-momentum 
tensor T^j^, in GR the Bardeen potentials satisfy 

= -AttGp, (5.3) 

= -47rG(p - 2V2 s) , (5.4) 

where Tqo = p is the energy density, and S the scalar part of the anisotropic stress tensor, 
see eqs. (A.9)-(A.I1). Therefore the conformal mode satisfies 

V2(cj + 47rGE) =0. (5.5) 

Imposing the boundary condition that the a and S vanish at infinity, the Laplacian is 
invertible, so we get a = —dyrGS. The conformal mode is therefore a constrained field, 
determined by the matter content, and vanishes if S = 0. 

The situation changes drastically in gravity. In this case, as we show in app. A, 
linearizing over flat space the conformal mode satisfies the equation 

(□ - ml) {a + 47rGS) = {p - 3P). (5.6) 

Thus the conformal mode a, shifted by 47rGE as in eq. (5.5), now satisfies a massive KG 
equation, where ttiq is just the same as eq. (2.25). Thus, the conformal mode becomes 
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dynamical, and describe the scalar particle that appears in Stelle theory. As we see from 
eq. (5.6), its source is given by the trace of the energy-momentum tensor, T = = 

— {p — 3P), as we expect indeed for the conformal mode. 

Consider now the theory with action 


S 



^P-Pai(D)p| , 


(5.7) 


where, at some UV scale p, ai(n) = ai is a constant, while in the IR ai(n) = In 

this case in the UV the conformal mode is described by eq. (5.6). In the IR, in contrast, 
the model reduces to eq. (1.2). The dynamical content of this theory, in the scalar sector, 
has been studied in Sect. 3 of ref. [7], where is has been shown that <I> and 'h are non- 
radiative field that satisfy Poisson equations, as in GR. The same is therefore true for the 
conformal mode a = (<I> — 'I')/2. We therefore conclude that a RG flows in which ai(n) 
evolves from a constant in the UV toward ai(n) = in the IR corresponds to a 

flow in which the conformal mode, which in the UV is dynamical, and is just the spin-0 
mode of B? theory with mass given by eq. (2.25), becomes a constrained field as we flow 
toward the IR. 


6 Dynamical mass generation for the conformal mode in 
Einstein gravity 

The discussion of the previous sections suggests that, because of IR effects, a mass term 
could be dynamically generated for the conformal mode. The introduction of non-local op¬ 
erators of the type Rn~‘^R allows us just to describe such a mass term, in a diffeomorphism- 
invariant manner. 

This leads us to investigate more closely the dynamics of the conformal mode in Ein¬ 
stein gravity, without R? corrections, to see if a similar “condensation” of the conformal 
mode might take place. Indeed, there are several indications for the relevance of the 
conformal mode in GR in the IR [80,81]. In classical general relativity the conformal 
mode is a non-propagating degree of freedom. At the quantum level, when gravity is cou¬ 
pled to massless conformally-invariant fields, it however acquires a dynamics through the 
conformal anomaly. As stressed in [81] , an anomaly implies that the effect of quantum fluc¬ 
tuations can remain large at the longest length and time scales, so quantum fluctuations 
can affect the behavior of the classical theory even in the far IR. This is most evident in 
D = 2 space-time dimension. In this case, the trace anomaly takes the form (see e.g. [41]) 

(r^) = —p, (6.1) 

where N = Ns + Np is the total number of massless scalar and Dirac fermion fields. The 
trace anomaly can be derived from the non-local Polyakov action^® 

5’anom = -^ J d^X^/^R^R, (6.2) 

^®The coefficient N in front of 5anom becomes V —25 if one also takes into account the metric fluctuations 
themselves, beside the fluctuations due to matter fields. 
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whose variation gives the right-hand side of eq. (6.1), and which can therefore be added 
to the classical Einstein-Hilbert action, to provide an effective action which takes into 
account the quantum effect of the trace anomaly. In L) = 2 one can always write the 
metric in the form (3.8), without truncating the theory, and in terms of the conformal 
mode the Polyakov action becomes local, 

jY r _ _ _ 

^anom = / (fx^/^{-(jn\a + Ra), (6.3) 


where the bars denote quantities computed with the reference metric In H = 2 

the effect of the anomaly is especially important, since the Einstein-Hilbert term is a 
topological invariant and the dynamics is entirely contained in Sanom- Observe that, in 
D = 2, the propagator of the operator is logarithmic, G{x,x') oc log(x — x')^, and 
therefore grows with distance. 

The situation is conceptually similar in four dimensions. In this case the trace anomaly 
is given by (see e.g. [41,42,82] for pedagogical introductions) 


(TH) 


hC^ + b' 



+ h''nR, 


(6.4) 


where, as in eq. (2.1), is the square of the Weyl tensor, E the Gauss-Bonnet term, and 
the coefficients 6, 6', h" depends on the number of massless conformally-coupled scalars, 
massless fermions and massless vector fields. The term on the right-hand side of 
eq. (6.4) can be derived from the variation of a local action proportional to f d'^x^—g R?. 
The remaining combination is obtained from the variation of a non-local action 


= -\ j (^E - A 


-1 

4 


b' I E-^OR] -2bC^ 


(6.5) 


where 

A4 = + 2i?^"V^n, - ‘^Ra + l(V^i2)V^ . (6.6) 

Once again, this non-local action becomes local in terms of the conformal mode. Writing 
the metric as in eq. (3.8), using as a fiducial metric g^i, = (and choosing appropriately 
the local R^ counterterm) one finds (see e.g. [80,81]) 

‘S'anom = j (□cj)^ , (6.7) 

where 

+ 62iVy — 28) -I- Qgrav • (6-8) 

Here Ns,Nf and Ny are the number of massless conformally coupled scalars, massless 
Weyl fermions and massless vectors, respectively, while Qgrav is the contribution from the 
gravitons, which is not unambiguously known, and —28 is the contribution from the a 
field itself [80,83]. 

Of course, a difference from the two-dimensional case is that in H = 4 writing the 
metric in the form (3.8) is a truncation of the theory. However, we see from eq. (6.7) that 


20 





the propagator of the a field in momentum space is oc 1 /k^, which in coordinate space, in 
D = A, gives again a propagator growing logarithmically, as in Z) = 2, 

G{x, x') = log - x'f] , (6.9) 

where /r is an IR regulator. This means again that the fluctuations of the a field are 
large at long distances, and justifies a posteriori the approximation of keeping only the 
conformal mode when studying the quantum effects of gravity in the IR. The fact that 
the conformal mode dominates the quantum theory in the IR is also confirmed by the 
fact that it gives the most infrared divergent contribution to the propagator in de Sitter 
space [84]. Observe that the fourth-order operator A 4 in D = 4 is the analogue of the 
second-order operator □ m D = 2, since it is conformally invariant, A 4 = A 4 . 

As discussed in ref. [83], the large fluctuations in the two-point function of the cr field 
due to its logarithmic growth create a situation that looks very similar to what happens 
in the two-dimensional XY model, with the Berezinsky-Kosterlitz-Thouless (BKT) transi¬ 
tion. Let us recall the argument of ref. [83]. The four-dimensional theory for the conformal 
factor has solutions of the form 


cj(x) = O'In jx — xo| , (6.10) 

which are analogous to the vortices of the 2-dimensional theory. The euclidean action for 
a, evaluated on this configuration, takes the value [83] 

, ( 6 . 11 ) 

where L is the IR cutoff and the “lattice spacing” a is the UV cutoff. The parameter xq 
is a collective coordinate corresponding to the center of the solution. The number of ways 
of placing this solution in a volume V is proportional to V, so the corresponding entropy 
is 4 log L/a, and the free energy is 

F = (IqV - 4) !■> (0 ■ (6.12) 

If > 8 /q the free energy diverges as we remove the IR cutoff {L —)■ 00 ) or the UV 
cutoff (a —>■ 0), so these configurations are irrelevant. In contrast, if < 8/q they 
dominate the partition function, so at the critical value of Q we expect a phase transition. 
For the solution with q = 1, = 8. The situation is indeed completely analogous to 

the BKT phase transition in two dimensions, which is also triggered by the logarithmic 
growth of the vortex-vortex interaction. The role of the coupling constant to be tuned to 
get the phase transition is here played by Q, which depends on the number of massless 
conformally-coupled particles. If the unknown graviton contribution in eq. ( 6 . 8 ) is not too 
large, the theory with just a single massless photon is in the region where < Q/., and 
the non-trivial configurations dominate. 

In ref. [83] it has been suggested that < Q/. corresponds to a phase where the space- 
time geometry is no longer smooth. Actually, pushing further the analogy with the BKT 
transition can suggest a different scenario. Indeed, in the two-dimensional case, in the 
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regime where the vortices dominate, the correlation length becomes finite and a mass gap 
is generated. It is plausible to expect that the same happens in the four-dimensional case. 
We therefore expect the generation of a term oc cr^ in the effective action for u, whose 
coefficient is given by a dynamically-generated mass scale. Of course, this mass term 
must be constructed in such a way that it does not spoil the underlying diffeomorphism 
invariance, which leads us again to a term see eq. (3.11). 

Adding this term to the anomaly-induced effective action (6.7), to quadratic order we 
get an effective action for a given by 


0 ( 2 ) _ _ 




The corresponding euclidean propagator is 

Stt^ f d'^k 

G(U = 


—ikx 


(27r)^ 


(6.13) 


(6.14) 


where = Ihvr^A^j^/Q 


Carrying out the angular integrals, we get 


2 /■'' 


dy- 


(6.15) 


where Ji(y) is a Bessel function. The final integral can in principle be carried out in 
term of a (not very illuminating) Meijer’s G-function. However we already see from this 
expression, together with the expansion Ji(y) — y/2 for small y, that the presence of the 
term regularizes the integral near the integration limit y = 0. In the limit —)• 0 

we recover eq. (6.9), except that the artificial IR regulator y of eq. (6.9) is replaced by 
the quantity y that appears in eq. (6.14), which now is a physical quantity related to 
Arr. In the opposite limit >> 1, introducing the notation u = (y^x^)^/^ and using 
as integration variable z = y/u, we have 


Gix) = 


2 

2 


1 


dz 


u 


1 + z^ 


^ I r 

V TT u3/2 Jq 


Ji{uz) 

^3/2 


dz 


1 + z^ 


Stt 


cos uz — 


(6.16) 


where in the second line we have used the asymptotic expansion of the Bessel function.^' 
The integrand is an oscillating function, which oscillates faster and faster as u increases, 
so the Green’s function vanishes faster than ^ i.e faster than l/(y^x^)^/^. 

The physical picture that emerges from these considerations is that, when the dynamics 
of the conformal mode is described by the action (6.7), fluctuations grow larger and larger 
at long distances, as shown by the logarithmic growth of the propagator in eq. (6.9). This 
is a sign of an IR instability of the vacuum state. When a mass for the conformal mode 
is dynamically generated, the long-distance correlations grow only until the separation 
|x —x'l reaches a length-scale of order 1/y, and after that they decay faster than l/|x —x'|^. 
Therefore, there is no IR instability. The dynamical generation of a mass for the conformal 
mode is the mechanism by which the vacuum restructures itself, as a response to the IR 
instability triggered by the anomaly-induced effective action. 

^^Observe that, as w —>■ 00 , the regime uz fixed corresponds to 2 ; —^ 0, whose contribution to the integral 
is suppressed by the factor in the integrand. The dominant contribution therefore comes from the region 
uz —>■ 00 , where we can use the asymptotic expansion Ji{x) ~ (2/7ra:)^^^ cos(a: — 37r/4) valid for large x. 
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7 Conclusions 


In this paper we have explored the possibility of generating non-local terms, such as those 
given in eqs. (1.1) and (1.2), through strong IR effects in gravity. We have hrst of all 
considered the running of the coupling associated to the B? term in higher-derivative 
gravity. We have stressed that a peculiar aspect of the running of this coupling is that, 
even in the deep IR, its calculation requires a knowledge of the UV completion of GR. 
Indeed, in GR the corresponding beta function is not even gauge-independent, and has no 
physical meaning. Stelle theory provides an example of an extension of GR where, thanks 
to the improved UV behavior, the beta functions of the 0{E?) couplings are physical 
and gauge-independent. However, the result from Stelle theory cannot be taken literally 
either, because the theory is plagued by a ghost (and, furthermore, existing computations 
of the RG flow have been performed only in the regime E ^ rn-p\). In the absence of 
a concrete reliable result for the RG flow, we have simply explored the consequences of 
a scenario in which the coupling 02 associated to the E? term grows in the IR. In this 
case, similarly to QGD, we expect the dynamical generation of a mass scale, Arr, and of 
a non-trivial form factor 02 ( 0 ). We have provided arguments, based on the analogy with 
the known non-perturbative behavior of the gluon propagator in QGD, that suggest that 
in the strong coupling regime the emergence of a pole-like behavior 02 ( 0 ) oc is quite 
possible. 

The emergence of a dynamical scale Arr would have far reaching consequences for 
solving long-standing puzzles concerning dark energy. In particular, this scale emerges 
from the logarithmic running of a dimensionless coupling constant, similarly to AqcD; 
and does not receive large loop corrections, contrary to what happens to the vacuum 
energy. It is therefore technically natural. We found that, to explain dark energy today, 
the dynamically-generated scale Arr should be of the order of a milli-eV. Of course we 
cannot predict its numerical value, just as we cannot predict the value of Aqcd. However, 
it is interesting that, to obtain a dynamical explanation of dark energy, we do not need to 
introduce a fundamental particle with an astonishingly small mass m Ho r\j 10 eV, 
nor a cosmological constant A with the amazingly small value A ~ Hq. In our scenario, in 
the IR the action is given by eq. (4.1), with a fundamental mass scale Arr ~ meV, which is 
by no means a particularly surprising value. The connection with the present value of the 
Hubble parameters comes out because, in a cosmological context, derivatives of the metric 
are of order if, so □ and □ ^ ^ H while R ~ Thus the term A^j^RD 

becomes competitive with m^^R when H ^rr/"^pi which, for Arr meV, gives H r-^ 
Hq ~ 10“^^ eV. However, in this scenario the fundamental scale emerging dynamically in 
the theory is Arr, and there is no fundamental particle with the astonishingly small mass 

10-33 gY 

Another interesting aspect of the above analysis is that it draws attention to the 
dynamics of the conformal mode of the metric, since the running of the coupling 02 as¬ 
sociated to R^, from a constant value in the UV toward 02 ( 1 ^) oc in the IR, has an 
interesting interpretation in terms of the conformal mode, as we discussed in sect. 5. This 
stimulates to look in more detail into the dynamics of the conformal mode in Einstein 
gravity, in which the higher-derivative term is provided by the anomaly-induced effective 
action, i.e. by the loops of the massless matter fields in the theory. Once again, the 
conformal mode seem to play an interesting role. In Einstein gravity supplemented by 
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the anomaly-induced effective action the emergence of a term R'0~‘^R can be interpreted 
as the dynamical generation a mass term for the conformal mode of the metric (an in¬ 
terpretation that becomes blurred in Stelle theory, where a mass term for the conformal 
mode is already present at tree level because of the R? term). We have argued that such a 
dynamical mass generation can take place as a consequence of the strong IR fluctuations 
generated by the two-point function of the conformal mode, which grows logarithmically 
at large distances and, as already suggested in in ref. [83], can trigger a phase transition 
similar to the Berezinsky-Kosterlitz-Thouless transition in two dimensions. If this results 
in a dynamical mass generation for the conformal mode, one again obtains in the IR a 
non-local model of the type (1.1) or (1.2). 
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A Dynamics of gravity in the scalar sector 

In this appendix we give more details on the computation leading to eq. (5.6). In GR the 
gauge-invariant scalar and vector perturbations of the metric are non-dynamical, and only 
tensor perturbations are dynamical. In R? theory this is different since, as we mentioned in 
Sect. 2.3.3, there is a scalar particle with a mass mo given by eq. (2.25). To understand how 
this takes place, and which scalar perturbation becomes dynamical, we linearize the metric 
around flat space, using the scalar-vector-tensor decomposition of the perturbations. This 
is the standard tool of cosmological perturbation theory, but can also be very useful for 
separating the perturbations over Minkowski space into pure gauge degrees of freedom, 
gauge-invariant but non-radiative components, and radiative degrees of freedom. We first 
recall how the formalism works in Einstein-Hilbert gravity, following [85,86], and we will 
then extend it to R? gravity. 


A.l Bardeen variables and conformal factor in GR 


Expanding the metric over flat space, the metric perturbation can 

be decomposed as 


^00 

hoi 

hij 


2V', 

+ da 

-2(t)6ij + (didj 


A -|- -{diCj + djei) + , 


(A.l) 

(A.2) 

(A.3) 


where /3* and e* are transverse vectors, dj/3* = = 0, and is transverse-traceless, 

= 0 and = 0. We have therefore decomposed the 10 components of into 

four functions {V', 0,7, A} which are scalar under spatial rotations, two transverse vector 
/3* and e*, which carry two independent components each, and a transverse-traceless tensor 
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hjj^, which also has two degrees of freedom^® The linearized massless theory is invariant 
under the gauge transformations —)■ It is useful to decompose 

also the gauge functions in the form = Bi + diC where = 0. In terms 


of these variables the gauge transformation reads 

(j)^(j) + {l/3)V^C, (A.4) 

\^X-2C, (A.5) 

Pi ^ f^i — Bi, ei ^ ei — 2Bi, (^-6) 

while is gauge invariant. One can then check that the combinations 

^ (A.7) 

T = i/)-j+^X, (A.8) 


are gauge invariant. Indeed, these are just the Bardeen variables, normally introduced 
when studying perturbations over a Friedmann-Robertson-Walker background, restricting 
to a scale factor a(t) = 1. Similarly, in the vector sector the transverse vector Hj = 
Pi — (l/2)ej is gauge invariant. Thus, out of the 10 components of the metric perturbation 
four are pure gauge degrees of freedom, corresponding to the two scalar functions 
A,C and the two components of the transverse vector Bi. The remaining six degrees 
of freedom are gauge invariant, and are organized into the two scalars and the two 
components of the transverse vector Hj and the two components of the transverse traceless 
tensor ■ Of course, not all these six physical degrees of freedom are radiative. To show 
this, we decompose also the energy-momentum tensor as 


Too 

= P, 

(A.9) 

Toi 

= Si + d,S, 

(A.IO) 

T- ■ 

= PSij + (didj — ^ T j + djUi) -\- aij , 

(A.ll) 


where cr® and 5* are transverse vectors and aij is transverse and traceless. Physically, 
p is the energy density, P the pressure, and S, crj and aij are, respectively, the scalar, 
vector and tensor part of the anisotropic stress tensor. These components are not in¬ 
dependent, because of energy-momentum conservation. To linear order over flat space 
energy-momentum conservation reads = 0. In particular, in the scalar sector, the 

equation = 0 gives 

V‘^S = p, (A.12) 

while the longitudinal part of = 0 gives [86] 

S = P+\v‘^T.. (A.13) 

o 

®®In the context of the scalar-vector-tensor decomposition, by scalar we will always mean scalar under 
spatial rotations. These quantities are of course not Lorentz scalars (for instance, ip is the (00) component 
of the Lorentz tensor hf^A- Similarly for the use of the words vector and tensor. 
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At the quadratic level in the action (i.e. at the linearized level in the equations of motion) 
the scalar, vector and tensor sectors decouple. Let us focus on the scalar sector. Choosing 
the gauge A = 7 = 0, we have = —cj), ip = 'i’, and the metric in the scalar sector can be 
written as 

ds^ = -{1 - 2^)dt^ + {1 + . (A.14) 

The quadratic part of the Einstein-Hilbert action, restricted to the scalar sector, is [ 86 ] 

= m|i j d^x (^-362 + . (A. 15) 

To linear order the interaction term is 

Sint = \j d^xh^,T>^\ (A. 16) 

which, restricting again to the scalar sector, reads 

= j d^x(3d>P + Tp) . (A.17) 


We can now find the equations of motion for the scalar variables, by taking the variation 
of the action with respect to <1> and The variation with respect 

to $ gives 


mpi 6$-2V2(<k-T) 


= -3P, 


(A.18) 


while the variation with respect to T gives the Poisson equation 


2mpiV2d> = —p. 


(A.19) 


Equation (A.19) shows that d* is a non-dynamical field, which vanishes if p = 0 (since, 
imposing the boundary condition that <1> vanishes at infinity, the inversion of the Laplacian 
is unique, and = 0 implies d> = 0). Equivalently, the initial conditions for $ are not 
free degrees of freedom, but are fixed by eq. (A.19) in terms of p, i.e. <b is a constrained 
variables. Since <b is uniquely determined by the matter source, the same is true for $ 
in eq. (A.18). Indeed, using eq. (A.19) we also have 2mpPV^^ = —p which, combined 
with eq. (A. 12) and with the invertibility of the Laplacian, gives 2mpj<b = —S. Using also 
eq. (A. 13), we finally get 

which, inserted into eq. (A.18), gives a Poisson equation also for T. Writing the result in 
terms of Newton’s constant, with l/m-pj = SttG, we get the standard result 

V2$ = -47rGp, (A. 21 ) 

= -47rG(p - 2V2s) , (A.22) 

^^Observe that, in principle, we should rather take the variation with respect to the original scalar metric 
component and 7 . However, because of the definitions (A.7) and (A. 8 ), the variation with respect 

to (j> is the same as that with respect to 4>, and the variation with respect to tp is the same as that with 
respect to I', while the variations with respect to A and 7 gives equations that are consequence of the first 
two, because of gauge invariance, see [ 86 ] for details. 
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Proceeding similarly in the vector and tensor sectors one finds 

= -IGirGSi, (A.23) 

= —WirGaij . (A.24) 

This shows that tensor perturbations are indeed radiative degrees of freedom, while the 
gauge invariant scalar and vector perturbations are non-radiative degrees of freedom that 
satisfy a Poisson equation, rather than a Klein-Gordon equation. In particular, in the 
absence of matter, = T = 0 and H* = 0 . 

It is also instructive to observe that, if we incorrectly truncate the theory in the scalar 
sector setting T = 0 in the action (A. 15), $ would seem a dynamical field, with a fully 
dynamical equation of motion obtained setting T = 0 in eq. (A.18), and furthermore its 
kinetic term in the action (A. 15) is such that it would seem a ghost! However, when 
one correctly keeps both T and <I>, one finds that it is the variation with respect to T, 
eq. (A. 19), that reduces 4) to a non-dynamical field. 

Finally, let us recast these results in terms of the conformal factor a and another 
independent perturbation ip, defined writing the metric perturbation, in the scalar sector, 
in the form 

ds'^ = [—(1 — 2ip)dt^ + (1 + 2 (^)dx^] . (A.25) 

At the level of linearized theory, = 1 + 2cr + Then, comparing with eq. (A.14) 

2a = - T , (A.26) 

2(p = <I> + T . (A.27) 

Since 4> and T are constrained non-dynamical fields, the same is true for a and ip. In 
particular, from eqs. (A. 21 ) and (A. 22 ), 

VV = -47rG(p - V^S), (A.28) 

VV = -dvrGV^S. (A.29) 

The second equation, together with the invertibility of the Laplacian assured by vanishing 
boundary conditions at infinity, implies a = —dvrGS. So the conformal factor is propor¬ 
tional to the scalar part of the anisotropic stress tensor, and vanishes in the absence of 
matter. 


A.2 Bardeen variables and conformal factor in gravity 

We consider now the dynamics in the scalar sector for the action 


■S'eh + Sji2 



1 

2 


rupiR — a2R^ 


(A.30) 


Since we are particularly interested in the dynamics of the conformal mode, and the latter 
drops from the Weyl tensor, which is conformally invariant, we do not consider the term 
here. Restricting again to the scalar sector and writing again the metric as in eq. (A. 14), 
to linear order in the perturbations we have 


R = 2 (^3$ - 2V24>-k , 


(A.31) 
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Then, the action to quadratic order in the perturbations, in the scalar sector, is given by 

c(2,s) _ c-(2,s) , c(2,s) ^{s) / . 

~ '^R'2 ^ ‘^int ) 

where “^int given in eqs. (A. 15) and (A. 17), respectively, and 


= -4a2 


d'^x 34> - 2V^$ + 


(A.33) 


We can now repeat the same analysis performed above for GR. The variation of the total 
action with respect to ‘h gives 

mh - 2V2(4> - T)1 - 802 ( 35 ^ - 2V2)(3$ - 2V^<^ + V^T) = -3P, (A.34) 


while the variation with respect to T gives 

2 mpi 4 > - 802(34' - 2 V^ 4 > + V^^') = - . 

Combining these equations we find 

mpi(4> + T) = S - V^p, 


(A.35) 

(A.36) 


and 


$ _ - 


m 


PI 


4> = 


1 


(p_V 2 s) + 


1 


-V- 2 p. 


(A.37) 


1202 J 3mpj • / ' 24a2 

Equation (A.36) is the same as the equation obtained in GR for the combination 4> + T, 
as we see adding eqs. (A.21) and (A.22), while eq. (A.37) is just a Klein-Gordon (KG) 
equation (with some source term) for a particle with squared mass mg = —mpj/( 12 o 2 ), in 
perfect agreement with eq. (2.25). Thus, the field ‘h becomes dynamical. Observe that in 
the limit 02 — 7-0 (at fixed mpi) we only retain the singular terms in eq. (A.37), and we 
recover the GR result, eq. (A.21).^° 

We can now rewrite the results in terms of the conformal factor a and of the variable 
(p, using eqs. (A.26) and (A.27). The equation for (p is again (A.28),^^ so this remains a 
constrained variable and satishes the same equation as in GR, while the equation for a 
becomes 


(□ — mg) (fj + dvrGS) = —— (p — 3P) 


(A.38) 


Thus the conformal mode a, shifted by dvrGS as in eq. (A.29), now satisfies a massive KG 
equation with mass mg. Its source, not unexpectedly, is given by the trace of the energy- 
momentum tensor, T = = ~{p ~ 3P). Thus, the scalar particle that appears 

in Stelle theory is described by the conformal factor of the metric. For 02 > 0 we have 
mg < 0, and in this case the conformal factor is tachyonic. 


^°This is completely analogous to what happens in linearized massive gravity, where again "h satisfies a 
massive KG equation [86-88] . It is interesting to observe that $ and are scalar only under rotation, but 
not under boosts, where they transform in a way dictated, respectively, by the transformation properties of 
hoo and htj, so they mix with the gauge-invariant variable in the vector sector. Similarly, even the action 
in the scalar sector is not invariant under Lorentz transformations, since boosts mix it with the action in 
the vector sector, which in turn mixes also with the tensor sector. The full set of transformations can be 
found in [86]. Nevertheless once the field $ becomes dynamical, its equation of motion nicely collapse into 
a KG form, as required by the fact that, in the end, its quanta becomes particles that must satisfy the 
relativistic dispersion relation. 

^^The equation for p will in general get an extra contribution from the term, which is also responsible 
for generating a dynamical field in the spin-2 sector. 
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